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Abstract 

We discuss possible tunneling phenomena associated with complex wave vectors along directions 
where the spin degeneracy is lifted in non-centrosymetric semiconductors. We show that the result 
drastically depends on the direction. In the [110] direction, no solution can be calculated in the 
usual way assuming that the wave function and its derivative are continuous. Conditions that 
permit one to find physical solutions are discussed and consequences are drawn. As a result, there 
is no spin filtering in such a direction but the spin undergoes a precession through the barrier with 
the rotation angle being proportional to the barrier thickness. In a direction close to [001] we find 
a spin-filter effect in close agreement with the model discussed by Perel et al. [Phys. Rev. B 67, 
201304(R) (2003).]. 
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I. INTRODUCTION 



Understanding spin-dependent tunneling through semiconductor barriers is a fundamen- 
tal problem in semiconductor physics. A description of this coherent process is crucial for 
spin-subband engineering of semiconductor heterostructures and superlattices. Moreover, 
spin-dependent tunneling through crystalline barriers has also become a topic of major in- 
terest in spintronics.^'^ This article thus lies at the interface between general semiconductor 
physics and spintronics. Regarding semiconductor physics, it has close connections with open 
problems in the envelope function theory.^ Harrison studied the problem of heterogeneous 
materials and introduced the conditions of discontinuity of the envelope function, taking a 
general viewpoint, well beyond the semiconductor area.- A decisive step in semiconductors 
was performed by BenDaniel and Duke^ who defined specific discontinuity conditions of the 
derivative of the envelope function between two media with different effective masses, based 
on the conservation of the probability current. This approach has been successfully applied 
to heterostructures by Bastard^ and has become the standard calculation routine, yielding 
very accurate energy positions of the energy bands, in perfect agreement with the experi- 
mental data.- Hereafter, analogously, we deal with periodic lattices which are perturbed by 
a spin-orbit potential and where, due to the absence of space-inversion symmetry, the spin 
degeneracy of the bands is lifted through a wave-vector-dependent "exchange" field.-'- We 
show that the matching conditions of the derivative of envelope function at the boundaries 
cannot be "as usual" . Thus, the basic tunneling equations are not known. 

First of all, dealing with tunneling phenomena requires an accurate knowledge of the 
energy structure in the forbidden band gap i.e. of the complex band structure of the bar- 
rier material. In pioneering articles, Heinel^' and Jones^i derived general properties of the 
evanescent states and showed their complexity over 6-dimensional wave-vector space, con- 
sisting of complex vectors associating a pure imaginary component to a real - propagating 
- one. It is often thought that the electrons will tunnel through such complex-wave-vector 
states as they would do through usual evanescent states (with pure imaginary wave vectors), 
and this intuitive explanation is probably supported by our familiarity with the tunneling 
of electrons located in semiconductor side valleys (e.g. in the conduction band of silicon). 
Hereafter, we deal with spin-dependent tunneling of conduction electrons through a gallium 
arsenide barrier, a compound with no inversion symmetry.^ Such processes were investi- 
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gated by Perel' et al. in a stimulating article, using the effective mass approximation and 
under simplifying assumptions; Quite large spin-filter effects were predicted.— Although the 
complex band structure of GaAs was expected to be well known, we have recently found 
that, in fact, the spin-orbit interaction and the absence of inversion symmetry had never 
been taken into account simultaneously throughout the Brillouin zone.— The evanescent 
band structure was calculated by several authors. Chang^^ considered semiconductors ori- 
ented in the [100], [111], and [110] directions, with space-inversion centers {Oh group) or 
without space- inversion centers (T^ group), but without taking into account the spin-orbit 
coupling. Chang and Schulman^^ performed a detailed calculation of the band structure 
of silicon, which belongs to the Oh group. Schuurmans and t'Hoofir^ studied semiconduc- 
tors belonging to the group but explicitly discarded terms which lead to odd k terms 
so that essentially they studied GaAs and AlAs as if they belonged to the Oh group. In 
Ref. Q the evanescent band structure in the fundamental gap of GaAs-like III-V semi- 
conductors, including both the spin-orbit coupling and the lack of inversion symmetry, was 
carefully calculated within a 14 x 14 and a 30 x 30 k ■ p Hamiltonian framework. Then 
it was demonstrated that the evanescent states in the fundamental gap present an original 
topology, with loops connecting opposite spin states at the center of the Brillouin zone.— 
This very structure has strong consequences for electron tunneling. Here, in order to remove 
any unnecessary complexity, we start dealing with electrons with a unique effective mass m, 
inside and outside the barrier - an approximation used in numerical applications in Ref. isl . 
The spin splitting in the barrier is described via the D'yakonov-Perel' (DP) Hamiltonian.- 
We thus revisit a classic in elementary quantum mechanics - the tunneling of free electrons 
through a square potential barrier - but in a case where the evanescent states in the barrier 
are spin split. From general considerations, we derive relevant boundary conditions which 
are sensitive on the crystallographic direction. We demonstrate that the tunneling process 
can become rather involved: The case of loop-shaped real-energy lines correspond to wave 
vectors which have both an imaginary component, which defines the tunneling direction, 
and an orthogonal real component so that one has to deal, so to say, with a "classical" 
tunnel effect in the sense where it is possible to recover almost usual tunneling properties 
- analogous to off-normal tunneling of free electrons - but in a subtle way. In the case of 
one-dimensional tunneling with a complex (neither real nor purely imaginary) wave vector, 
the tunnel effect seems to be "anomalous": A spin precession occurs around a "complex 
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magnetic field". We show that the derivative of the envelope function, which is the solu- 
tion of the Schrodinger equation, undergoes discontinuities at the barrier plane - usually, 
in semiconductor heterostructures, discontinuities of the derivative arise as a consequence 
of the different effective masses in the well and in the barrier material^ - and we propose 
a treatment of heterostructures. After entangling the two spin channels, it is possible to 
recover a situation which has strong analogy with standard tunneling and where the discon- 
tinuity of a "magnetic current" can be viewed as the result of a kinetic-momentum transfer 
at the barrier interfaces. The spin-orbit-split barrier exerts a torque on the electron spin, 
similar to spin-torque phenomena in ferromagnetic junctions as predicted by Slonczewski^ 
and Berger— , but in this latter the barrier is constituted of magnetic material, a 

spin transfer occurs between the tunneling electrons and the magnetization. 

The layout of this paper is as follows: In Sec. [Tll we give the background relative to 
the spin splitting and to the conservation of the probability current which will be used 
afterwards. We show how the spin splitting can lead to complex (not strictly imaginary) 
wave vectors in the barrier and we analyze the consequences on the probability current. In 
Sec. Ill, we study a barrier normal to [110], and in Sec. IIVI we look in detail at the case of 
an incident wave whose direction is almost normal to a [001] barrier. A summary is given 
in Sec. |Vl 

II. BACKGROUND 
A. Symmetry 

Let us present the notations used throughout the present article (see Fig. 1). e is a unit 
vector. The direction of the axes, defined by e^, e^, with respect to crystal axes, will be 
given in each case, is normal to the barrier, = eno in Sec. Ill and = eooi in Sec. IIVI 
We define k/= ^ + q, kjj = + Q+iK (^, q, Q, K are all real vectors), and kiii= ^ + q. 
We also introduce the following notations: p =xex + yey, ^ =C,^ex + C,yey, q = qe^, Q =Qez, 
K =Kez, k/-r = k///-r =^^x + + qz, kjj ■ r =^^x + ^yU + (Q + iK) z. Without spin, the 
wave function of the incident plane wave and in the barrier should be written as e**-^ ''"''^^-' 
and e^[i-p+'yQ+^i<)A respectively. 

To describe the structure of the evanescent states, we use the k • p method. In a n- 
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band model, the energy dispersion curves result from the diagonalization of a (ri x n) k ■ p 
Hamiltonian H, but k is a complex vector so that H is no longer Hermitian and the evanes- 
cent states are associated only to real eigenvalues To find the energy dispersion 
curves, we have to solve the secular equation detM(k) = det[H — EI], where I is the 
identity. Because the Hamiltonian is Hermitian when k is a real vector, we have the re- 
lation M(k)** = M(k*). Thus, detM(k*) = [detM(k)*]* = [detM(k)]*. It follows that 
EnQi) = -E„'(k*), where the band indices n and n' may or may not refer to the same 
band.— >ii Moreover, Kramers conjugates correspond to the same energy, so that the state 
associated to (k, \up)) and the state associated to (— k, \down)) are degenerate.— Let us 
recall that Kramers-conjugate states are obtained by application of K, the time reversal 
operator: K = —iayKQ where ay is the relevant Pauli matrix and Kq is the operation of 
taking the complex conjugate.— Thus in GaAs, the spin degeneracy is lifted and we expect 
that the four states [(k, |s)) , (k*, \s')) , (— k*, |— s)),and (— k, |— s'))] be degenerate, |s) 
and \s') being up-spin states in directions which, generally, are not parallel (Fig. 2). We 
are going to see a concrete example in Subsec. IllBl where \s) and \s') are quantized in the 
same direction, and in Sec. HV] where \s) and \s') are not quantized in the same direction. 

B. Energy levels 

In the introduction, we mentioned that the evanescent band structure is deeply altered 
when the lack of inversion symmetry is taken into account together with the spin-orbit 
splitting. A particular topology consisting of loops connecting Kramers-conjugate spin states 
near the zone center was shown along directions of the type K[^/K, 0, i] when the ratio 
^jK = tan 6* is fixed. Such loop structure can be expected to arise as it is known that 
a band cannot stop.^° Depending on 6', we obtain the different pictures shown in Fig. 3. 
The spin vector along a loop is defined by the mean value of the Pauli operator ct. In the 
small-fc limit, we get two opposite spin vectors. When going off the zone center, a numerical 
calculation shows that the two spin vectors rotate to become parallel at the point where the 
two subbands are connecting. The appearance of these loops is the fingerprint of a strong 
band mixing of the first conduction band and of the three upper valence bands with remote 
bands (more precisely with the second conduction band). Indeed, as long as the wave vector 
remains in some vicinity of F, the energy levels are well described by the DP Hamiltonian, 
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where the spin states in the subbands only depend on 6 (see Sec. HV]) . Observe, in Fig. 3, 
that - because the extension of the loop tends to zero when 6 tends to 45° - the portion 
of the loops which can be described in this analytical model also has an extension which 
can become vanishingly small. Hereafter, we stay in the framework of the DP model, which 
allows analytical calculations. 

Throughout the present paper, we take the origin of the energy at the bottom of the 
conduction band so that the relevant Hamiltonian is written as 

H = Hq + Hdp 
Hdp = 7 X • ^ 

where m is the effective mass. Hdp is the DP Hamiltonian which describes the spin- 
splitting:-^ = (k) = [x^, Xy, Xz] = [kx {kl - kf) , ky {kl - kl) , k^ {kl - fc^)] . When k 
is real, the energy levels are pure spin states, quantized along in the plane perpendicular 
to k. Note that the two eigenvalues of % ■ are opposite, equal to the square roots of 

= + + Xz- We designate by x+ (X-) square roots of x^ (x^with a positive real 
part and x_with a negative real part, if relevant). x+ (X-) "^iH be used in Eqs. 12.61 and [2751 
The eigenvalues of H are written as S (k). 

Inside a finite-width barrier, the incident plane wave e*''^ is usually to be replaced by 
e^^^ which corresponds to an imaginary wave vector -^iK. 

a) If the incident wave vector k/ is in the [001] direction (k/ = [0, 0, g]), the wave vector 
in the barrier is k/j = [0, 0, ±.iK\ and the degenerate eigenvalues of H are 8 (k) = — 7^/^^ 
which is the (real) energy E (k) in the forbidden band gap. If k/ is almost in the [001] 
direction (kj = [^, 0, q] with ^ ^ g), k// = [^, 0, ±.iK] and the eigenvalues of H are 8 (k) = 
— 7j, (i^^ — ^^^K\J — ^'^ which is the energy E (k) in the forbidden band gap as well. 

b) If k/ is in the [110] direction (hj = ^ [110]j, a simple idea would be to take k// = 
±^ [110] which leads to 8 (k) = —-y^K"^ ± i-^K^. This quantity is not real and cannot 
be an energy E{k.)M We are therefore led to consider a wave vector such that k// = 
^(Q±zi^) [110]. 

The calculation is given in Appendix A. The resulting band is plotted in Fig. 4, over a very 
broad energy domain to reveal its general structure. We are only interested in evanescent 
states located in the forbidden band gap, i.e. states with a small negative energy. For our 
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purposes, a key point is that, at a given energy, we have exactly the four possible states, the 
wave vectors of which are {Q±iK) with spin | and {—Q±iK) with spin |, the latter being 
obtained from the former through K. In short: 

(k) = (k*) = E^ i-k) = E^ i-k*) (2.2) 

Eq. 12.21 provides us with a concrete example of the ideas developed by Jones^i who 
showed that E (k) = E (k*). The corresponding four plane waves are e'^^^'^^^^ |, e'^i-Q^^^) | 
or e^^^e*'^^ |, e^^^e^**^^ | (this is schematically shown in Fig. 5). This leads us to define 

^= e'^' T ^= e-'^' i (2.3) 

so that the four plane waves write e^^^ ft and e^^^ 4- 

In the following, | and | are the up and down spins when the x vector, which plays the 
role of a magnetic field, lies along a real direction and is taken as quantization axis. When 
X is not coUinear to any real direction, the spin eigenstates are tk and l^: In Subsec. IIV Aj 
we shall see that |k and |k are no longer orthogonal. The implications of a wave vector 
k = Q ± iK in the [110] direction will be considered in detail in Sec. III. 

C. Probability current 

1. The free-electron probability current 

We consider a spin-orbit-split barrier separating two regions where the electron states 
are described by plane waves and where the potential is taken equal to zero, as shown in 
Fig. 1. The barrier potential is assumed to be a positive constant. When dealing with 
tunneling phenomena through crystalline barriers, the wave vector component $, parallel to 
the barrier plane has to be conserved. For an incident plane wave, which has a real wave 
vector component parallel to the surface plane, this implies that the imaginary component 
of the wave vector inside the barrier has to be orthogonal to the barrier plane. Then, the 
imaginary component of the wave vector inside the barrier defines the tunneling direction. 
To analyze the tunneling processes, we distinguish two different mechanisms: i) The wave 
vector has collinear real and imaginary components along the normal to the barrier (we 
refer to this mechanism as ^^ara-type) and ii) The real and imaginary components of the 



7 



wave vector in the barrier are orthogonal (we refer to this mechanism as ortho-type). We 
would point out that a plane wave with the real (imaginary) wave vector (^ + Q) (resp. zK) 
is associated with the "classical" probability current, e.g. calculated for a free electron, 
3-^ = ^(^ + Q)/m (resp. 0). Such currents, with a zero divergence, conserve the local 
probability in any domain located in the barrier. On the contrary, a plane wave with the 
wave vector (^ + Q) + iK is associated to J-^ = e~^^'^ h {$, + Ql) / m : It looks as if the 
local probability were to be no longer conserved in a domain located in the barrier, unless 
Q = 0, because V • J-^ = — {2h/m) (K • Q) e^^^'^ . The case Q = 0, results in a laminar 
free-electron probability flux. The loops in the complex band structure which have been 



studied in Ref. 
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correspond to ortho-tunneling, the normal to the barrier plane lying along 
[001], a direction where the spin splitting is zero. On the contrary, tunneling along the [110] 
direction, a direction where the DP field is maximum, is a para-process. More precisely, the 
definition of the free-electron current probability, 

m = Re 



m 



h 

— Im[V^*VV'] (2.4) 
m 



is obtained from the conservation of the local probability when the potential in the 
Schrodinger equation is real.— Obviously, the equations expressing the conservation of the 
probability have to be re-examined carefully in our case, where the Hermitian potential is 
non-real due to the spin-orbit interaction. The detailed derivation of the relevant current 
operator which allows one to calculate the true currents of probability J-t and J = J+ + J„ is 
given in Appendix [Bl There, it is shown how to extend the usual procedure, which consists 
to define the velocity v from the relation 

^ dH , , 

2. The ortho- and para-processes 

Coming back to the specific case of the GaAs-type barrier, let us derive a few basic results 
and introduce some definitions. The orbital part of the wave function of the conduction 
band is S in usual Kane's notation^ and we write ip^ = S ] (k) and = 5* | (k) where 
t (k) =t [see i) below] or T (k) =tk [see ii) below] and i (k) =i [i)] or i (k) =ik [ii)]. The 
corresponding Schrodinger equation is: 



i) Ortho-process. Let us assume that k/j = ^+iK ( i.e. Rek//. Imk// = 0, with ^ ■ K = 0) 

is a possible evanescent state. Because E = (^^ — K"^) is real on a real-energy line, the 

2m 

terms x±(k) originating from the spin part of the Hamiltonian are also to be real. We follow 
the usual procedure to derive the expression of the probability current, tk and |k are no 
longer orthogonal but in any case the real spin term disappears so that we obtain 

|2 

L = _v . J/ J (2.7) 



dt 



which is the usual relation for probability conservation. Care has to be taken that the 
relation V ■ 3^ [^±1 = V ■ J [4'±\ does not mean that J-^ [4'±\ = J [V'i] • However, in such 
a case, a number of classical results derived for free electrons will be recovered. 

ii) Para-process. In the case of one-dimensional tunneling along the n direction, where n 
is a unit vector normal to the barrier, which involves a complex wave vector k = (Q + iK) n, 
x(k) = {Q + iK)^ X (n)) we quantize the spin along the direction of x ip) which is a real 
non-normalized vector. x±(k) are no longer real. We follow the same procedure to derive 
the expression of the probability current and we obtain 

I 1 2 

d\i'+\ r r 1 2 I i2 , . 

^^ = _V.j/[^J+-^7lmx± IV^il (2-8) 

These equations could suggest an interpretation in terms of two-channel transport with 
a generation-recombination rate, analogous to Giant MagnetoResistance phenomena.— In 
such a case, we would classically expect a spin mixing and we will show that, indeed, a 
formal analogy exists. However, care has to be taken that, at a given k, ipj^ and do not 
correspond to the same energy except when x is zero. 



3. The [110] direction 

More specifically, we will deal with electron tunneling along the [110] direction, a direction 
where the spin splitting is maximum in the real conduction band. On this example, we 
illustrate the preceding considerations. Let us consider for instance the up-spin channel, 
where a possible wave vector is k = (Q + iK) eno as shown in Sec. Ill Bl with the wave 
function 

^P^{z) = e^(«+^^)^ (2.9) 
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and the DP field 



The free-electron current is 



x+ = - (g + iKY 



3+ \'ip+\ = — Im'?/;1VV'+ = — 



-2Kz 



On a real-energy line (see Eq. IA5I1 



2h 



m 



KQe 



-2Kz 



(2.10) 

(2.11) 
(2.12) 



dt 



2K 



1 



'2Kz _ g 



(2.13) 



Along the real-energy line, the eigenstates of the Schrodinger equation comply, as ex- 
pected, with the continuity equation, the current J+ being to be identified. Here, it is easy 
to show that (see Appendix [Bl) 



^ / 



(2.14) 



In the real conduction band, taking ip^ = e^'^^, we obtain 



27c ,37 2 



1 5 



2 I 

7c? ± o79 



hdq 



E{q) 



(2.15) 



2h hdq 

Concerning an evanescent wave = e'^^^'^Q)^ ^ it is easy to check that 3± \Qi^^'^Q)^^ = Q 
on a real-energy line. 



4- Waves conserving the free-electron probability current 

The waves which conserve the free-electron current of probability play a special role: They 
appear to be "quasi-classical states" which allow us to build solutions yielding intuitive 
physical interpretations. The waves involved in an ortho-process verify Rek ■ Imk = 
and we have seen in Sec. Ill C 21 that this condition ensures the conservation of 3^ . In 
the case of a para-process, a paradigm being tunneling along [110], J-^ is not conserved in 
a given spin channel. Therefore, it is necessary to consider an intricated wave function, 
ip (r) = ip+ (i") T +4'- (r) i= V^4- T +4'- i= T i- III the following, we indifferently 
use the notation ip-^ and ■?/'_ or and ip^ 
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The free-electron probability current is given by^° J'^ [V^] = (1/''^) R-s ('0*p'?/')o. 
where the index a means a summation (partial trace) on the spin or J-^ = 
(l/m)Re [4'\.p4'_^ + 4'*_pi'_) = J-^ [i^+l + J'^ b^-]- -^^^ Kramers symmetry, the wave 
functions in the barrier '?/'//+ ^^^^ i^ii- can be written 



The free-electron probability current carried by the function of the type 
(Aae-^^ + 526^^) e'^^^ is (e = ±1) 

J/ [0] = A [2KluiAlB2 + tQ (2ReA25* + |A2|'e-2^" + |52|'e2^^)] 

In the barrier, let us write = ^1^//+ t +^//- i) so we have 

h 



(2.16) 



(2.17) 



[^//]=-<!2Q 



+Q 



+ 2K 

\A2 



lm52A2 + Im52v42 

2' 



A, 



Bo 



\B, 



We see that the probability current in the barrier is constant if and only if 1^2 1 

B, 



(2.18) 



A, 



This leads to A2 = Ae'^^, A2 = Ae''^^, B2 = Be'^^, and B2 = Be'^^^ where A 
and B are two complex numbers. So the general expression of a wave sustaining a constant 
3-^ inside the barrier is 



It is useful to write 



(2.19) 



(2.20) 



where 

Sx = Sx{z)=Xii +X* ^ (2.21) 

The Kramers conjugate of S\ is S\ = KSx where K is the time reversal transformation. 
Observe that S\ and Sx are eigenstates of the helicity operator (p ■ ct) for the eigenvalue 
hQ. 

Let us look at the spin direction defined by S\. Recall that the spin quantization direction 
is along the x (^iio) vector. We call Oz' the direction parallel to x (^iio); Ox' and Oy' are 
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in the 11-;^ plane normal to x (eiio)- The spin direction is defined via (cr^/), ((7^/), (cr^/). First 
of all we note that {uz') — while {ux') = 2ReA^ and {ay>) — — 2ImA^ for S\ (0). The spin 
is in the 11^ plane. Any spin direction in the 11^ plane, that we call an in-plane direction, 
can be described by a suited value of A. For instance with A = ex^iOx^ {o'x') = cos26\^ and 
{(Tyi) = — sin26'A, apart a common factor, 9x being the angle between the Ox' axis and the 
spin direction. 

It can be shown that the largest vectorial space consisting of J •''-conserving waves at a 
given energy is 6 = {^a,/?} where 

^ = (a^e-^^ + pBe''') t + (a^^e"^" + /3*Be^") ^ (2.22) 

with a and /? e C. ^ is a vectorial space over M, but not over C. 

Moreover, the existence of a superposition principle implies that any linear combina- 
tion with real coefficients of two solutions with a current of probability of a given sign 
has to be a solution associated to a current of probabihty of the same sign. This is 
a strong constraint which is verified over ^o= {^a,b} ® {Sa} = {{Ae'^'^ + Be^^) Sa], 
a vectorial subspacc of € (in this subspacc J-^ [^A,BSa] = 2 |q;|^ J-^ [$_4^g]), or in 
^{^a,bS)x — cos 9 ^^^j3 Sx -I- sin6' ■§^^a,b {^^>^ | - at fixed 9 - which also is vectorial 
subspace (in this subspace 3^ {^a,bS)x e = 2 |q;|^ J-'' cos 2^). 

D. Standard tunneling case 

The standard tunneling case is to be recovered when 7 is zero, therefore, we build our 
analysis in close relation with it. A crucial point is that the probability current has to be 
constant so that R + T — 1 where R {T) is the refiection (transmission) coefficient. 

We shall need the standard (without spin) function ■0'-°^ (z) defined as: 



V'^"^ {z) 





(^) 


= 4 


- hie^^"' 


(z<0) 








+ 626^" 


(0 < 2; < a) 




(^) 






(a < z) 



(2.23) 



where ^ < 0, < z < a, and a < z respectively correspond to the incident wave (index /), to 
the wave in the barrier (index //), and to the transmitted wave (index ///), as illustrated 
in Fig. 1. 
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tp^ (z), a C^-function, meets the boundary conditions 

,(0). N ,(0)^ ^ ^^"(2:0-) dlp^^'^ {zo+) 

i/j"- > (2:0-) = i)^ ' (zo+) , = -Q-^ , 2:0 = or a (2.24) 



hi 2 (g2 + K'^) sinh Ka exp A'a»i (g^ + K"^) 



ai D (q + iKf 

a2 2q{q + iK)e^°- expKa->i ^ q 

ai~ D ~ (g + iK) 

62 _ 2g (-g + iK) e'^" expXa»i ^ g(-g + ii^) _2;^^ 

03 ^ exp£a»i . giTe-'^" j,^ 

ai D ^ ^\q + ^K)^'' 

D = {q + iKf e^'^ - (g - ii^)' e"^" 



(2.25) 



The function ijj^'^'' [z) is such that the probabihty current J-^ 



(0) 



is constant. The 



reflection coefficient R = \bi/aif and the transmission coefficient T = 1 03/01!^ are such that 
R + T = 1. 

Also observe that, if we multiply ip^^^ by any C^-function /(r, |, |), the new function 
and its derivative are continuous at the interfaces, satisfying the initial boundary conditions. 
Consider the case where the incident wave is e^'^'^. If we take /(r, f, |) = e*^'"" t or /(r, f 
, i) = e*^"" I, we obtain a solution to the tunneling problem if, and only if, the incident 
component e*^'*"'"^^ '" and the reflected component e*'-"'*"'"^-' '", correspond to the same energy.— 



III. PARA PROCESS: [110]-ORIENTED BARRIER UNDER NORMAL INCI- 
DENCE 

A. General considerations 

In the case where the wave vector is parallel to the [110] direction and the Hjjp Hamil- 
tonian is taken into account, we have seen in Sec. Ill BI that the wave vector is to be of the 
form {eQ ± iK) eno to get a real eigenvalue (an energy) of the Hamiltonian. But in such 
case, J-^ is not conserved (Eq. I2.17P and even not constant inside the barrier, so that the 
standard calculation routine to flnd the solution (i.e. the continuity of the wave function 
and of its derivative) cannot apply. 
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We could try to build a solution according to the usual procedure, but with a wave in 
the barrier involving the two spin channels, which can give a constant 3^ (see Eq. I2.19P 



^i{z) 
^11 iz) 



(3.1) 



where fl' and |L are defined in Eq. 12. 3[ 

The usual boundary conditions {C^ function) for the down-spin channel for instance yield 
four equations determining Bi, A2, B2, and A^ 



qBi 

I2 {Q - iK) e-^W-^-^)« + ^2 (Q + iK) e-^(Q+*^)'^ 



B2 + A2 

{Q - iK) A2 + {Q + iK) B2 
A^e'^" 

-A^qe'^" (3.2) 



They only provide a non trivial solution if the determinant of the system is equal to 
which gives the relation 



(g2 -Q^ - K^) sinh Ka + 2iKq cosh Ka = 
The only solution is K = but it is not relevant to our problem. 



(3.3) 



B. Solutions to the tunneling problem 

1. The constant-^ case 

We go back to the Schrodinger equation to determine the proper boundary conditions and, 
to avoid any unnecessary mathematical complexity, we here assume that 7 is constant over 
the three regions. Along the [110] direction, with k = (l/y^ A;[110], the DP Hamiltonian 
writes 

Hdp = J,k' ± ^jk^ (3.4) 

where the + (resp. — ) sign applies to the up (down) spin, quantized along the DP field. As 

d 

usual, we obtain the effective Hamiltonian by substituting k with — iV , i.e. k with —i-^r-- 

az 

Hdp = -7.^ ± (3.5) 
14 



Thus, we have the two equations 

-7c 
-7c 



92 


1 93 1 








1 93 " 







= \E-V{z)]i^^ 



(3.6) 



where ^ (z) = when < 2; < a and (-2) =0 outside. Because the DP Hamiltonian 
was obtained using the perturbation theory, we will look for a solution to the effective 
Schrodinger equation to the first order in 7 only. Let us consider the up-spin channel: We 
write 

= + (3.7) 

where ip^^^ is the standard function (obtained for 7 = 0, and defined by Eq. I2.23j It is a 
function, with a discontinuous second derivative). V^j^^ is a first-order term in 7 so that the 
Schrodinger equation to the first order writes 



7c 



+ 0^7" 



[E-V{z)]^^ 



" dz^ ' 2^' dz^ 
We integrate this equation from one side of the interface to the other 

-I zo+e 1 r / ;^2,/,(0) ' 



(3.8) 



7c 



d%jj 



dz 



ZQ-e 



1 

+ 2' 



7 



dz"^ 



zo+e zo+e 



[E-V{z)]^^dz 



(3.9) 



ZQ-e 



ZQ-e 



Then 



lim i -7c 

e — >0 



dz 



-\ zo+e 



ZQ-e 



1 

+ 2' 



7 



zo+e • 



zo~e. 







(3.10) 



Taking the standard function (Eq. 12.231 ) and referring to the limit at Zq inside the barrier 
and inside the well respectively cLS Zq and z^ , we obtain 



outside the barrier, and 



dz"^ 



dz"^ 



^0 



(3.11) 



(3.12) 



inside the barrier. At the interfaces i/j^^^Zq ) = ip 



^(°)(2o), then 



-Z7 



dz"^ 



^0 



^0 



(3.13) 
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This provides us with the jump of the derivative at the interfaces. To the first order in 
q/K 



dz 



^0 



1.7 



^ Ic ^ Ic 



(3.14) 



after Eq. 

Similarly, for a down spin Qi = —Q^, and we have 



dz 



^0 



^0 



(3.15) 



It is worth remarking that this very discontinuity condition was found in a quite different 
situation, involving Rashba-split quantum wells.— 

Now let us assume that = Q. The wave function constructed from the eigenstates in 
the three regions is 



iIjj(z) = Aie''^' + Bi{q,K,Q)e~''^' {z < 0) 

il)jj{z) = A2{q,K,Q)e~^'e'^' + B2{q,K,Q)e^^e'^' (0 < ^ < a) (3.16) 

i,,jj{z) = A,{q,K,Q)e^'^^ {a<z) 



with the coefficients Bi{q, K, Q), A2{q, K, Q), B2{q, K, Q), and A^^q, K, Q) to be determined. 
To the first order in Q, the solution can be expanded as 



-iqz 



(z) = {a,e"^' + b,e-"^') + (3,Qe 
^jj{z) = (ase-^^ + 626^") e'«^ + Q (ase^^^ + /^se^") e'«^ (3.17) 



with 



a2 



dA2{q,K,Q) 



dQ 



dB,{q,K,Q) 
dQ 

; P2 = 



and 



"3 



dAs{q,K,Q) 
dQ 



dB2{q,K,Q) 
dQ 



Q=o 



We write 



-0 = + V2' 



(3.18) 
(3.19) 

(3.20) 

(3.21) 
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where 



(^f {z) = aie^«^ + bie-'i^ {z < 0) 

(p^jj (z) = 036^5^ e'^" (a < 



(3.22) 



ipl (z) = /^iQe-^?^ < 0) 

^ ) ^s^ ^ g (ase-^^ + /^se^^) e^'?^ (0 < z < a) (3.23) 

v 

ip^ is a continuous function but its derivative is not. To the first order, its jump at the 
interfaces is 



dz 



(3.24) 



^0 



As we have derived that the jump of the derivative of the wave function ip is 2iQil)fj {zq), 
we deduce that (p^ is a continuous function and that the jump of its derivative at the 
interfaces is 



dz 



(3.25) 



^0 



This provides us with the four equations which determine the four coefficients P^, 02, 



and a'n 



iqpi - Ka2 + Kp^ 





#ff(0) 



(3.26) 



The solution of this system is 



Pi — —— 036*^" sinh ai^T = —ai sinh a/T 
K D 



-,Ka 



OL2 



a2 



P^ = ioge^^"- 



2K q + iK 



2K q-iK 



«3 = 



(3.27) 



In the following, we consider Q/K — ^K/A'j^ but also q/K as first order terms and we 
look for solutions up to the first order. The term in the reflected wave function arising from 
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Qfil is a second-order contribution which has to be neglected. Note that, in Region /, if the 
incident wave has the wave vector g, the reflected wave should have the wave vector — g', 
where q = qo — Sq and q' = qo + Sq. From Eq. IB27t it can be verified that 6q = (7/47^) q^ 
= [qo/KYQ. Then 6q is a second-order term which has to be neglected so that media / 
and /// have no sizable spin splitting. This indicates that the solution we obtain in the 
case of a constant 7 also constitutes a plausible physical solution when 7 is a step function, 
with 7 = outside the barrier. Also note that, at this level of approximation, is a wave 
which only exists inside the barrier and is not coupled to the free-electron waves outside 
the barrier. Because A3 = 036**^", we see that there is a pure dephasing between the up- 
{Ql = Q) Ei'iid the down- {Q^ = —Q) spin channels. 

We have to be sure that, in our treatment, the probability current is conserved along 
the tunnel process. The wave in the barrier, in the up-spin channel, is of the form ip [z) = 
{A2e-^' + B2e^') e'^' = <p {z) e^^^ with A2 = 02 (l - i%) and B2 = 62 (l + • Let us 
calculate J to the first order in Q by making use of Eq. 12.141 



J [^] =J. 



+ 



1_ 

2h 



d_ 

dz 



dz'' 



(3.28) 



It is sufficient to evaluate the term in the bracket to the zeroth order, substituting i(j with 
One finds 

2 



J M 



dip 



(0) 



dz 



[ip] + 2 



hKK 



dip 



(0) 



dz 



dz"^ 



iP^'^^iP^'> 



iP 



(0) 



(3.29) 



[iP] - 2^Q 



iP 



(0) 



2 (02^2 + ^262 



with 



3' 



\ m dz I h \ dz I h 



iP 



(0) 



27c 



Im ip 



.dip 



(0) 



h \ ^ dz 
By comparing these expressions, one obtains 



27c I 



+ {alh2 + a2hl) + 

n h 



iP 



(0) 



(3.30) 
(3.31) 



n 



Im ip 



dz 



iP 



(0) 



(3.32) 
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This definitely establishes current conservation in the tunnel process. 
Starting with an incident state Iv^o); transmission asymmetry T in the spin-dependent 
tunneling process as can be expressed as 



l|T(|^o))ll'- 






2 


l|T(|^o))ll' + 






2 



In the present case, we find T = 0. Whatever the incident spin, the tunnel barrier acts 
as a pure spin rotator, without any spin filter effect. 

The cases of a spin-split quantum well confined between infinite walls and grown along 
the [110] direction is discussed in Appendix B2. 



2. A unified description 



Let us now consider transport in the real conduction band, in Region / or ///. In the 
case 7 = 0, the solution of the Schrodinger equation is ip (z) = ip^^^ (z) = aje'^'^"^ + 6je~*'^°^, 
j = 1 or 3, 63 = 0. When 7 is non zero, the wave function, in the up-spin channel, has to 
be of the form 



%jj {z) = e 



i-d 



aAl + a—] e*««" + ( 1 + /3— ) e"*^"" 



-iSqz 



(3.34) 



where e*'^ is a phase factor. Here again, let us calculate J to the first order in 6q by making 
use of Eq. 13.281 Substituting ijj with in the bracket, one obtains 

2 



2h 



h 



dip 



(0) 



dz 



(0) 



- 2 (a*6,e2^««" + a,h*e-^"'^') 



iP 



(0) 



+ 4^5g(|a,f Rea- |6,f Re/5) - 2^5q 



>7c 



(0) 



(3.35) 
(3.36) 



For Re a = — Re /3, one finds 



Jf 



Jf 



+ 4^5gRea(|a,f + - 2^5q 

n n 



,7c 



+ 4^5gRea 
n 



iP 



(0) 



iP 



(0) 



(3.37) 
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Taking Rea = — 1 and Ima = lm.(3 = 0, 

J [ifj] = 

In the barrier, we consider (cf Sec. Ill B 1) 



(0) 



(3.38) 



a2 I 1 - 1 e""^^ + &2 ( 1 + I e 



K 



■Q ] ^K. 

K 



.Q 



iQzQ 



dz 

In the well, let us take (cf Eq. 13.341) 



ipw (z) = e 



iQzg 



72 ^V'g^ (^0) 
472c 9z 



-i5qj{z-zu) 



(3.39) 

(3.40) 
(3.41) 



(3.42) 



where zq is the boundary relevant to the region of the well (e.g. = in Region / and 
= o in Region /// ). Although we are still dealing with a unique effective mass and a 
constant 7, for the subsequent discussion, it is convenient to refer to 7 (7^) as 72 (72c) or 
7j (7^-^) where j = 1 or 3 in the different regions, and to the wave vectors as gjo — ^Ij and 
-{qjo + ^qj)- 



(zo) = e'Q^^ 
We obtain 



iQzo 



,2 



^B(^o)-V^iy(^o) = e^^^° 



JQzQ 



, 72 ^^g^ (^0) , 7, c}^!^^ (^o) 
472c i^^; 47^-^ 9z 



(^0) 



(3.43) 
(3.44) 



.,(0) 



JQzo 



^iQzo 



i 


5^2^ (^0) 


( ^2ljc 


^•72c 


47,c 


dz 


V 72c 


7jc 


i 






7i72c 


472c 


dz 


V 72c 


7jc 



(3.45) 



Here, we have used the relation 



(0) 



72< 



dz 



7 



dz 



(3.46) 
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which originates from the usual relation expressing current conservation in the absence of 
DP field.- When 7 and 7^ (i.e. m) are constant, ip^ {zq) — ip-^ {zq) = 0, which establishes 
the continuity of the wave function. 

Now, let us examine the matching conditions of the derivative 



K 



K 



(3.47) 



dz 



',(0) 



dz 



+ 2^QV'^°^ (^0 



(3.48) 



dz 



g^g^ (^0) 

dz 



- 2^5^,^^°^ (^0) 



6q, 



z5g,.^(°) (^0) 



(3.49) 



72c 



dipsizQ) dipw{zo) 



dz 



dz 



d^f (^o) djj'^ (zq) 

72c ~ 7jc 



,(0) 



dz 



dz 



(3.50) 



This is exactly the jump of the derivative calculated in Eq. 13.141 up to the second-order 
terms. Thus, starting from the standard solution, we have constructed in a very simple 
way a wave function which is continuous, associated to the constant current of probability 



(0) 



and which is the solution to the tunneling problem. 



3. Insight into the step-function case 

The case where 7 (z) = 'jg {z) is not a constant raises difficult questions. The problem is 
not to solve Eqs. 13.61 but to define a proper Hamiltonian, which has to be Hermitian: This 
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would not be the case simply by substituting 7 with 7 (z) in these equations and there are 
several ways to symmetrize this Hamiltonian. This is in line with the BDD approach when 
dealing with an heterostructure where m = m{z), i.e. where m depends on z, for instance 
m = mi in Region / and m = 1712 in Region 11.^'^ In that case, the starting point is the 
Hamiltonian 

The key idea is to transform this equation by defining the BED Hamiltonian 

where v is defined in Eq. 12.51 Then, an integration of the Schrodinger equation around the 
origin, exactly as performed above, will allow us to show that J = iIj*viI! is continuous because 
ip and ^^^^ are continuous. The BDD Hamiltonian guarantees probability-current conser- 
vation and the problem receives sound foundations. Unfortunately, the more complicated 

2 

term - 



form of the current of probability given in Eq. 12.141 - in particular due to the 
makes an analogous transformation not obvious so that the general case still remains an open 
question. However, let us point out that, when the masses and the DP-field coefficients are 
not very different over the three regions - a frequent situation in heterostructures - through 
the procedure described in the preceding subsection, we are able to construct a wave which 
is continuous at the boundaries and that conserves the current of probability. Therefore, 
this wave is a plausible solution. The principle is first to solve the envelope-function problem 
in the absence of DP field, i.e. g (z) = 0, taking into account the mass discontinuities in 
framework of the BDD formalism. This determines the standard function ip^^^ (z). Sec- 
ondly, the wave functions in the different regions are modified according to the rules defined 
in the preceding subsection (Eqs. 13.391 and 13.421) . The current of probability remains equal 



to 



(0) 



in the three regions. Concerning the continuity of the wave function at the 
boundaries, we have at zq = 

di^^^ (0) 



dz 



2iqiai 

5q ^ ^-^^ql and Q ^ (3.53) 

727ic 472c 472c 



Thus (see Eq. ^M\f 



^5 (0) - (0) = — ^ ^ ) (3.54) 

^ 72c \ 72c Ijc 
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We introduce 



^3,2 = and 57^-2 = 

7,- = 7,-,2 + <^7,-,2 72 = 7,-,2 - <^7,-,2 



(3.55) 



r.-,2 



7jc + 72, 



_ Ijc - lie 

2 ' --h-^- 2 



and 5r,2 



77C = Tj, 2 + ^Fj- 2 72c = rj,2 - S^ir 



(3.56) 



(0) - i^w (0) 



1 gitti 



2 72c 



1 - 



(5ri,2 



1 - 



<5ri,2 

ri,2 



-72 



7i,2 ri,2 
i2c "V 7i,2 ri,2 

At 2^0 = a, the situation is similar with 



1 - 



1 + 

ri,2 



, 91 Q (hl,2 JTl,2 



(3.67) 



dz 



(3.58) 



In the case where and are small and considered as first-order terms, the discon- 



^3,2 



tinuities are third-order terms which can be safely neglected. 



4- Quasi- classical picture (Regions I and III without sizable spin splitting) 



In the case where Regions I and /// have no sizable spin splitting, wc develop a quasi- 
classical picture of the tunneling process. For an up spin Q — Q^, so that the wave function 
in the barrier writes as 



ii+y 



K 



K 



iQz 



(3.59) 



The wave function for the down spin is obtained by replacing Q by —Q. We can combine 
the two spin channels to build the quasi-classical solution (z) corresponding to an incident 
wave with a spin lying in the plane perpendicular to the DP field 



ri (z) = (A T +A* i) (z) = s.m^r (z) 



',(0) 



(3.60) 
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which yields 



Defining 



i^ni^) = i^fh^)Sx{z) - ^^§-/ni^)^KS,{z) (3.61) 



tane = %^e (3.62) 
K 



we can write to the first order 



ruiz) = cos^<;(z)5a(z) - siTi9^^^^pf]{z)iKS,^ {z) (3.63) 

The transmitted wave is 

V;^,, = {Xe'Q^ T +A*e-*«» i) ase^"^ = Sx{a)^fn (3-64) 

The incident wave corresponds to a spin lying in the plane, normal to x(eiio). An 
important result is that the transmitted wave has the spin iSq, (a), i.e. rotated by the angle 
-2Qa. We can estimate the angle 2Qa ^ 0.2 [k/1 A j (a/1 A) in GaAs along the [110] 
direction, the largest reasonable value of K being smaller than 0.1 A~^, a value beyond 
which the spin-splitting in is no longer valid. The spin-split harrier appears to exert a 
spin torque which produces a rotation of the spin of the transmitted electron around the 
quantization axis, which is the direction of the DP field. There is no spin transmission 
asymmetry . The spin-orbit-split barrier acts as a spin rotator inside the 11^ plane. This 
has some analogy with the refiection of a neutron beam on a ferromagnetic mirror discussed 
in Ref. I28I which physically results from spin precession during the time spent by the 
evanescent wave inside the barrier. But, in this example, this straightforwardly arises from 
the difference in the reflection and transmission coefficients for the two spin eigenstates. 
Anyway, this spin precession provides an estimation of the tunnel time r, by using this 
built-in Larmor clock.— The effective fleld is determined through: Ml 27 |x| whereas 
VLt = 2Qa ^ {aiXe/lcl We flnd r ^ \ah/2'j^\ \Xe/x\ K"^ ■ In the [110] direction, Xe = 1/2 
(see Sec. [nB]) so that r ^ \ah/4:-f^K\ ^ IQ-^^ (a/1 A) (lA~^/K^ s. 

We recognize that the in-plane solution belongs to the subspace of free-electron-current 
conserving waves studied in Sec. Ill C4[ In that sense, we have restored a "classical" tun- 
neling process. Note that J = J-^ is a constant, but the classical "magnetic current" in 
Region II , 5J-^(z) = J^{z) — jfiz) is not and undergoes a discontinuity at the boundaries. 
Quite generally for any two-component spinor ip = ip, | | with ip, = $ e**^^ and 
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h 



h 



53f = -lm 

m 



T i) W T i) = - Re [^^ (p ■ a) 



m 



Thus, the jump of ^J-^ is 



Im 



m 



More exphcitly, we have 63 j = 63jjj = 0, 63^jj{z) ~ 2 |A| 
viewed as a kinetic-momentum transfer along the internal-field direction during the tunnel 
process, in strong analogy with the spin transfer resulting from spin torque in ferromagnetic 
structures, as introduced by Slonczewski and Berger.-ii^. 



m 



Re [^t (p.^)^];o 



r/ 



r/ 



m 



1 d^niz) 



K dz 



(3.65a) 
(3.65b) 
(3.65c) 



(3.66) 
. This can be 



IV. ORTHO PROCESS: [001]-ORIENTED BARRIER UNDER ALMOST NOR- 
MAL INCIDENCE 

It is not possible to stay in simple band schemes - like in Fig. 3 - as ^ has to be conserved: 
The relevant scheme is drawn in Fig. 6. To simplify without altering the physics of interest, 
the component of the wave vector normal to [001] is taken parallel to [100]. The spin is 
quantized along the Oz axis, taken parallel to [001]. As shown below, the eigenstates of the 
spin are in a direction normal to Oz. The energy writes 



E = -7,(/C2 - e) ± l^JC^JC^ - e (4.1) 

[E + lA^' - e)] ' = {iiK,f {le - e) (4.2) 

where the generic wave vector is ^eioo + ^/Ceooi- This equation may admit 4 real roots, ±i^, 
and -^K' . The states of the four wave vectors (,^, 0, zi^) tk, (^,0,— z-ft') tk*, {i.^.^K') |k', 
(^, 0,— i-fr') Ik'* have the same energies: K and K' are such that E ] (K) = E I {K'). 
Note that Kramers conjugate states, which would involve — ^, are not relevant because ^ is 
conserved. We introduce Kq = {K' + K)/2 and 6K = K' ~ K (note that this definition 
differs by a factor of 2 of the definition used in Sec. Ill, where 2Sq = q' — q The choice made 



in the present section makes more easy the comparison with the results derived in Ref. 
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13, the incident- wave 



We assume that K' > K > so that 5K > 0. Moreover, as in Ref. 
energy is smaller than half of the barrier energy, which means that q < K. 

As recognized by Perel' et al., the tunneling problem admits simple solutions under the 
approximation ^/Kq <^ 1. Besides, the spin asymmetry which originates from the spin-orbit 
interaction is characterized by the ratio 6K/Kq, which, from band-structure calculations^^ 
and from spin-precession experiments,—*^ is known to be small, i.e. 5K/Kq ^ 1. We 
further assume that aK^ is not small compared to unity, which corresponds to a barrier of 
small transparency, and consequently we have exp (—20 7^0) ^ 1. These three quantities, 
C,/Ko, 6K/Kq, and exp {—2aKo) will be hereafter taken as first-order quantities and we will 
look for solutions to the first order only. This does not imply that the quantity a6K = 
{aKo){6K/KQ), which is of crucial interest as it characterizes the spin selectivity of the 
barrier (as illustrated by the simple evaluation indicated below), is smaller than unity. In 
the physical problem, we consider electron tunneling under off-normal incidence and the 
angle of incidence is significant only when q and are of the same order, which means 
q/Ko <^ 1. We shall use this additional approximation only when it will be necessary to get 
analytical expressions of the wave vectors (Subsec. llVCj) . Intuitively, if we start with an 
unpolarized electron beam, the up- (down-) spin electrons merge from the barrier with an 
amplitude of probability almost proportional to exp —a K (respectively exp —aK') so that 
the current asymmetry - which, in this case, is also the polarization 11 of the current - is 
given by 



^-2aK _ ^-2aK' 



^-^n? ^ = t&nhaSK (4.3) 



Indeed, in Ref. 
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it is found that the polarization V of the transmitted current, when 
the primary beam is not polarized, is P ~ tanha^i^ (see below Eq. I4.15p . In practical 
cases, a6K cannot be larger than a (often small) fraction of unity. Nevertheless, in the 
calculation, we do not put any restrictive assumption on a bK (which is not assumed to be a 
first-order quantity) and we will calculate eigenvectors, when required, as a power expansion 
in a 5K, but, obviously, we keep in mind that the first-order term will generally be sufficient 
to reach a reasonable accuracy. 



26 



A. Zeroth-order wave functions 



The wave vectors K and K' are related through the equation 
{K' > K and assuming 7 > for the sake of simphcity) 



7,(^2 - e) - liKjK^ -e = -i.{K" - e) + i^k'Jk'^ - e 



(4.4) 



or 



7,(i^'' - K^) = 7^ (^i^V^' - + - 
Up to the first order in 6K/Kq, Eq. 14.51 writes as 



(4.5) 



(4.6) 



or 



We now calculate the eigenvectors. Let us write 

k = (^, 0, rjiJC) with r] = ±1, IC = K 01 K' , ^, K and K' >0. x = ^^{K, 0, ir]^). 

'X.z 'X.x ^•^y 



(4.7) 



The eigenvalues of a ■ x 



are ±^1C^K? - 



To the first order in E,/Kq, the normalized eigenvectors Ci | +C2 



are such that 



wave 
vector 



































iK 








iK' 







(4.8) 



spm X 



1 + 



2K 
2K. 



2K' 



2K' 



1 + 



2K' 



2K' 



Observe that |k and |k are not orthogonal (even in a first-order calculation - compare 
the first term to the third one after substituting K' with K). Inside the barrier the wave 
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function is of the shape \l/^^(r) = e''^'f "^^^ (z) and 



2K 



2K} 

+ A2 



1 + 



1 - 



1 - 



1 + 



2K 
2K\ 



2K' 
-(1 + ^) 



1 + 
-(1- 



2K' 
2K' 



;(-i2e-^'^ + 52e^" 



+ 



1 

1 

K'z 



Outside the barrier, we are looking for the solution of the shape: 





1 




1 




1 




1 




1 




1 


e"*^" + Ai 


-1 




-1 



-iqz 



and 





1 




1 




^3 












1 




-1 





The wave function writes as 



^!"{z) 



^>'"{z) 



[Aie"'' + Bit-"''] [1 1]* + 



[1 - 1]* 



A,e-^^ + 526^^ + ^ (-i2e-^'^ + ^26^^'^ 



\l 11 



+ 



2/s: 



(^26-^^ - 526^") + i2e-^'" + 826 



K'z 



[1 



[Age'^^j [1 1]* + A^e"^' [1 - 1 



(4.9) 



(4.10) 



(4.11) 



(4.12a) 

1]* (4.12b) 
(4.12c) 



The continuity of the wave function (Eqs. I4.12p and of its derivative at z = and z = a 
provides a linear system of 8 equations. A full discussion is given in Appendix O . This 
calculation has strong similarities with Slonczewski's^ approach of the tunneling between 
two ferromagnets separated by a barrier, because we deal with two coupled spin channels. 
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B. Polarization 



The transmission asymmetry T is 



T 



\t+f + \t-f 



2 I 4- 1 2 — 2 

with It"*" I = \A^\ (calculated when Ai = 1, Ai = 0) and \t~\ 



A-, 



(4.13) 



(calculated when 



Ai = 0, Ai = 1). All the coefficients Aj and Aj are calculated in Appendix [Cl 
To the zeroth order in ^/Kq, t"^ = t^, T = Tq. Now 



and we get the result of Ref. 



13 



iqKe' 



-Ka 



namely 



AqK'e 



l„-K'a 



{K' - iqY 



(4.14) 



Tn = tanh a 6K 



Up to the ffist order in ^/Kq , = tf, T = Ti, [tf 



^+ |2 



Ai 



but 



(4.15) 



and 



I are of second order in ^/Kq so that, up to the first order in $,/Kq, the result is the 
same as for the zeroth order: Tq = Ti. 

It is easy to show that this transmission asymmetry is nothing but the spin polarization 
of the transmitted beam when the primary beam is unpolarized, % = %_ = V. As we have 
only assumed that q < Kq, we may wonder why the ratio q/Ko does not appear in V. The 
answer is given if we perform the calculation one order further in 6K/Kq <^ 1. Then, a 
lengthy calculation leads to 



V 



tanha6K + 

Kp+q Kq 



(4.16) 



In the limit where SK/ Kq is negligible, V = tanh a 6K is recovered. 

Let us consider the transmission of a primary electron beam with an initial current po- 
larization Pi through a spin-filtering structure characterized by the transmission coefficients 
^-2aK' ^^-2aK^ (dowu) - spiu clcctrous. As the incident up- (down-) spin current is 

proportional to 1 + Vi (respectively 1 — Vi), the current polarization of the emerging beam 
is simply given by V 

[1 + Pi)e-2-^ - (1 - V^)e-^''''' Vi + Ii 



V 



-2aK 



+ (l-7^,)e-2»^' l + IiVi 



(4.17) 
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where 11 is given by Eq. 14.31 The above formula yielding the polarization of the transmitted 
beam is a standard expression for spin filters (in spin polarimetry, 11 is referred to as the 
Sherman function).— Thus, V in Eq. 14.161 appears to result from the combination of a 
primary-electron-beam polarization Vi ~ —5K/Kq when q/K^ <^ 1, which does not depend 
on the barrier thickness, with the spin asymmetry of the material, 11 = tanh (a (57^). The 
initial polarization —5K/Kq could be straightforwardly understood as resulting from the 
band mismatch, an interface effect. If this analogy provides us with a useful physical insight, 
it must, however, be realized that the above calculation is only valid when exp aK^ ^ 1 and 
cannot be extrapolated to a = 0. In any case, it is clear that Pi builds up in the early stage 
of the transport process. 



C. ^/Kq first-order wave function 

It is shown in Appendix ICl that there is no ^/Kq first-order term in ^2, A3, and -62- 
We are therefore going to calculate E,/Kq first-order terms in y42, -B2, and ^3. To be 
consistent with Subsec. IIV At we assume that 7^ and = 0. We obviously have to 
invert the role of K and K' if we start from Ai = and Ai 7^ 0. Let us recall that the 
calculation is performed with 5K/ Kq <^ 1 which is always true and S,/ Kq <^ 1. 

Eqs. Elil [Cni give 



(4.18a) 



-A2(l + ^le 



-K'a + 5^ ( 1 _ ^ , g 



11 e^'« = — 
K' / 2K 



K 



2K' \ K ' 



(4.18b) 



The determinant of the system defined by Eqs. 14.181 is 



Det 



iq 



-K'a 



1 , 



K'a 



(4.19) 



which differs from zero, therefore A2 and B2 can be calculated. 
We assume a 7^ (the case a = has no interest) and we obtain 



Ao 



2K' 



sinh (K'a) sinh (K'a) 



(4.20) 
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and 



Bo 



2K' 



sinh (Koa) 



(4.21) 



sinh (K'a) sinh (K'a) 

Noticing that i) ^/K' = ^/Kq (1 + 6K/2K) ^ {^/Kq) (1 - 6K/2Ko) ^ ^/Kq (the same 
result holds for ^/K ^ ^/Kq), ii) aSK < aK, iii) A2 cx Ai (Eq. ??), and iv) B2 oc 
Aiexp(-2ira) (Eq. ??), we get 



A2 
B2 



2Ko 



sinh (K'a) 



2K. 



Aoe 



-j^oa Sinh(Q'^^/2) ^ ^^^^g^j^/a sinh (iToa) 



sinh {K'a) 



sinh {K'a) 



(4.22a) 
(4.22b) 



From now on we assume that expi^oa ^ 1 so that sinhi^'oo/ sinhi^^'a = exp —a 6K/2 



and 



(4.23) 



A lengthy calculation shows that: 

i) Bi is proportional to {^/Kq) {5K/Kq) and therefore is negligible. However, we can 
note that Bi is not strictly equal to zero so that the reflected wave has a [1 — 1]* component 
even though the incident wave has only a [1 1]* component. 

ii) 



B, 



2K. 



iK + q _„i,Ki2 . . a^K 

2 -e "^-^/^smh hi 

iK-q 2 



B, 



We furthermore assume that q/K^ <^ 1 so that 



Bo 



2Ko 



-a&K 



[2 



-aSKl 



e B2 



and eventually 



A, 



li ( abK 
smh 



- 2 sinh^ 



abK 



A, 



(4.24) 



(4.25) 



(4.26) 



2Kq \ 2 
There is no assumption on abK in Eq. I4.26[ 

We note that, as A3 differs from zero, the incident wave with only a [1 1]* spin component 
is transmitted with a component along the [1 — 1]* spin direction. This means there is not 
a pure spin-filter effect along the x-quantization axis.— 



V. CONCLUSION 



Electron tunneling in a semiconductor with no inversion symmetry and in the presence 
of spin-orbit coupling involves complex wave vectors in the barrier. In directions where the 
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D'yakonov-Perel' (DP) field is non zero, the problem becomes highly non-trivial. We have 
distinguished two particular types of tunnel processes: Para-type process where we have 
one-dimensional tunneling with a complex wave vector and ortho-type process associated 
with a complex wave vector with orthogonal real and imaginary components. For a para- 
process, the DP field is a complex vector but it remains coUinear to a real direction so that 
the eigenvectors are orthogonal spin states. We have shown that, along the [110] direction 
no solution exists. The expression of the current of probability is re-examined, proper 
boundary conditions are derived, and a treatment of heterostructures is proposed. Quasi- 
classical states are shown to be in-plane solutions, which imply a pure spin rotation of the 
transmitted beam around the direction of the DP field. In the [110] direction, there is no 
spin-filter effect. This contrasts with the situation in the real conduction band where the 
spin-splitting is maximum along [110]. For an ortho-process, the DP field is a complex vector, 
which is not collinear to any real direction, and the eigenvectors of the Hamiltonian are no 
longer orthogonal spin states. Moreover, the evanescent eigenvectors are not associated with 
the same spin depending whether they propagate from left to right or from right to left. In 
this case, we have derived a first-order solution to the tunnel problem, which has strong 
similarities with standard off-normal tunneling, and an almost pure spin-filter effect was 
demonstrated, a conclusion consistent with the result of Perel' et al.— whose expression for 
the transmitted polarization has been corrected by the introduction of an initial interface 
polarization. 

All these questions should now be addressed experimentally and we think that experi- 
ments are within reach. For instance, further developments of the study of the polarization 
of a reflected spin-polarized electron beam can be considered, in line with the measurements 



reported in Ref. 



35l . Polarized-luminescence experiments in quantum wells grown along 



the [110] axis could also bring valuable information, as well as measurements on resonant- 
tunneling devices or photogavalnic-effect measurements in coupled quantum wells.— '^'^'^ 
The results derived in the present article provide new insight in spin-dependent tunneling 
in solids whereas they also open stimulating perspectives for spin manipulation in tunnel 
devices. 
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APPENDIX A: EVANESCENT BAND IN THE [110] DIRECTION 

Let us write k = (Q + iK)e, having in mind e along the [110] direction: e = eno = 
^ [110]. We have to find the relation between Q and K to get a real eigenvalue of the 
Hamiltonian H. This real eigenvalue is the energy. The Hamiltonian H writes as 

^ = 7c (Q + ^Kf + 7 ^-X = 7c {Q + ^Kf + 7Xe (Q + ^Kf ^ -e^ (Al) 

= x/ 11x11 (provided ^ 0). Xe > ^ dimensionless parameter, depends on the 
direction. If e = eno, X is parallel to e^o with Xe = 1/2. 
The eigenvalues are 

E (k) = 7c (Q + ^Kf + exe7 (Q + ^Kf (A2) 

The spin is quantized along 76^ so that 67 > corresponds to the spin ] and 67 < 
corresponds to the spin |. Separating the real and imaginary parts of the eigenvalue, we 
obtain 

Re^: (k) = 7c(Q' - K^) + %7(Q' - ^QK^) (A3) 
Im £ (k) = 27,gK + eXe7(3Q'i^ - K^) (A4) 
Looking for the real-energy lines, we have the equation 

lm£(k) = =^ 27,g + eXe7(3Q' - i^') = (A5) 

-2 Oyn2 I O ^ / , . ,1 7c 



= 3g' + 2e^Q = + e 4^Q if e = eno (A6) 

7Xe V 7 



Eq. IA6I is the relation between Q and K we were looking for. The energy is 

7Xe 



-e47Q^ - 87,Q2 _ ^A^q jf g = eno (AT) 

7 



For a given -E'((5) value, we have two possible choices of K, 



K = ±^3Q^ + 2e^Q (^=±^3g2 + e4^Qif e = enoj (A8) 

Let us note that |e4 (7/7J Q\ > SQ"^ so that \Q\ < \K\ and 

K ^ ±v/(4e7c/7)Q (A9) 
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The sign of £7 determines the sign of Q (7^ > 0). As stated above 67 > 0, which corre- 
sponds to spin I, gives Q > whereas 67 < 0, which corresponds to spin J,, gives Q < 0. 
We have the symmetry property 



E±{Q) = E^{-Q) (AlO) 

The study of the function E{Q) is straightforward and we take e = — 1 in the following, 
the other case being deduced by symmetry: 

d^-(Q) _ 24x,7Q' - 167cQ + {= I27Q' - l^lcQ + 4^ if e = eno^ (All) 



dg ^^'-^ 7Xe V " 7 

The roots Qi and (52 of the derivative are 



Incidentally we note that 



27Xe 67Xe 

Qi = -,Q2 = ^ife = eno (A12) 
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£;-(Qi) = (A13) 



The corresponding curve is plotted in Fig. 4. It must be realized that we are only dealing 
with evanescent states, which correspond to a negative energy. Thus, for a given energy 
£J < 0, we have two possible Q-values, ±(5, each associated with a given spin subband. 

Finally, we find that, at a given energy, we have exactly the four possible states, the wave 
vectors of which are {Q ± iK) with spin f and {—Q ± iK) with spin [, the latter being 
obtained from the former through K. In short 

(k) = (k*) = E^ (-k) = E^ (-k*) (A14) 

APPENDIX B: CONTINUITY EQUATION AND DEFINITION OF THE PROB- 
ABILITY CURRENT 

1. Definition of the probability current 

Consider a Hamiltonian given by 

H = ^ ttjPj + ^ bjkPjPk + CjkiPjPkPi + V (Bl) 

3 3,k j,k,l 
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where aj, bjk, Cjki are Hermitian matrices, invariant under permutation of the indices i, j, k 
, and where V is real. We define the velocity operator 



dH 
dpj 



% + 2 ^ bjkPk + 3 CjkiPkPi 



(B2) 



k,l 



It will be useful to take the following notations: 



^1 (r) T +^2 (r) i= 



V'l (r) 
V'2(r) 



(P, (r) T +02 (r) 1= 



02 (r) 



(r) (r)] , ((/)| = [0t (r) 0; (r)] , (^IV') = 01 (r) V'l (r) + 0^ (r) ^^2 (r) 

^1 (r) ^1 (r) + (r) (r) = 1^ (r) |' = |^|' 
P^i (r) 
pV'2 (r) 

[p*0l (r)] [pV'i (r)] + [r0; (r)] [pV'2 (r)] 

(r)] [p^i (r)] + [-p^l (r)] [pt^^ (r)] (B3) 



= [p*0I(r) ^0; (r)] = [-P0t (r) -p0;(r)] 



The Schrodinger equation is 



ih 
—ih 



dt 
dt 



^ ajpj IV') + &jfcPjPfe IV') + XI ^jkiPjPkPi IV') 



The continuity equation can be written 



"'•Is'") + (I*!* 



dt 



J2 [(V^lajPjV^) - (PjV^lajV^)] 
j 

+ Yl ii'^l^jkPjPk'iP) - {PoPkij\bjkij)] 



+ X [i^l^jklPjPkPl'^) - iPjPkPl'il^lbjk'^)] 
j,k,l 



Note that 



or 



{Pjij\ajtp) = {ip\ajPjipy 



(B4) 



(B5) 



(B6) 



(B7) 
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Similarly 



ipjPkM^jk'^) = {MhkPjPk'^T ; iPjPkPi'iPlbjk'iP) = {McjkiPjPkPii^T (bs) 



Therefore 



^ (^MOjPjV') + XI ii^hkPjM) + XI ii^\(^jklPjPkM) 
■ 3 j,k j,k,l 



(B9) 



The probabihty current J has to satisfy 



V • J = -- Im 

h 



X + X i'^l^jkPjPk'^) + X i'^\(^jklPjPkM) 

. 3 j,k j,k,l 



^ vjw + v-j(^) + v-j(=^) 



(BIO) 



From the expression of the velocity operator, we tentatively define the probability current 



as 



- {ip\ajip) + X ii^lbjkPki^) + 9 X (Pki^l^jkiPk'ip) 



k j,k,l 

where cc refers to the complex conjugate. We calculate 



+ cc 



3 



Let us consider the first term 



= ^ i^M) +cc^ ii^M) 



(Ell) 



(B12) 



(B13) 



X ^ l^^^ i'^M) = i X ('^\^3P3^) - 



3 3 



= --ImX(V'|%P,-V') = V-jW 

3 

The second term gives 

jf ^ = X (^l^JkPk-ijj) + CC = X [(V'l^ifePfc^) + (PfcV'l^ifcV')] 



(B14) 



(B15) 



X Vjjf^ = ~ iPMbjkPk'ip) + iPk'iPlbjkPj'ip) - iPjPk'iPlbjki^)] 



3,k 



2 

Im y ^ ii'ibjkpjpki') = V • J 



(2) 



(B16) 
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Concerning to the third term 



J(3) _ 3 



kd 



kd 



3i 



■ [(Pk^l^JklPjPl'^) - iPjPk^lCjklPl'^)] 7^ V • J^^^ 

k,l 

Let us now consider the quantity 

^PjPkPl (V'lCjfeiV') = Y [{'^\CjklPjPkPlll^) - {PjPkPMCjkl'^)] 



(BIT) 
(B18) 



jkl 



jkl 



-^^[{Pj'^lCjklPkPl'il^) - {P3Pki^\CjklPli^)\ 



jkl 



= YpjYpkPi{'4^\cjki'4^) 



kl 



h 



Y^j^PkPl (V'ICjfeiV') 



(B19) 



kl 



We have 



^PkPi (V'IcjwV') + J 



^(3) 

i 



k,l 



Im J] (^|c,«p,M^) = V • J(=^) (B20) 



j,k,l 



Thus, we can define 



r(3) _ t(3) 



^PkPl (^|CjW^) 



k,l 



Finally, the j component of the probability current can be taken as 



or 



- {ip\ajip) + ^ (iplbjkPki^) + 2 XI {Pk'^\cjkiPk'^) + 2 ^PkPi ii^lcjki'ip) 



= + + 3 ^ {Pki>\CjklPk'(p) + ^PkPl {ll^\Cjkli') 



(B21) 



+ CC (B22) 



(B23) 



kd 



2. Quantum well grown in the [110] direction 

To illustrate some simple consequences, wc apply the preceding results to the practical 
case of quantum wells grown in the [110] direction. First, let us point out that, in this case, 
a direct calculation of the current of probability is straightforward 

,dijjj^ d'^ip^ 1 9^V± 



dt 



= -7, 



dz'^ 2 dz^ 



± 



"9^2 ' 2 ' dz^ 
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(B24) 



Multiplying the first equation by 4'*±, the second equation by ip± and subtracting them, 
we obtain 



dt 



-Ic 



(B25) 



or 



d\ 



dt 
-V- 



dz^ 



d-z"^^ 



9\ ^ ,2 



(B26) 



We consider a well made of a spin-split semiconductor (GaAs) confined between infinite 
walls located at z = and z = a. At energy E, for a given spin, the wave function (p {z) 
consists of a combination of eigenstates associated to the wave vectors q (E) and —q' (E) 
(see Fig. 6, upper part) which satisfy 



/3 



(B27) 



The wave function writes 



= Ae^«" + 5e~*^^ (B28) 
and verifies the boundary condition: 0(0) = (j){a) = so that A = —B and q + q' = 



or 



A straightforward calculation gives 

|2 r, 



2tA sin (^^) e-^^«^ 



V J 



dt 



-- \Af sin {{q + g')^;} 



1 







(B29) 



(B30) 



due to the energy expression (Eq. IB27p . The probabihty current J is conserved as it should. 
However, a calculation of J according to Eq. IB26I yields 



a2 



(B31) 



Obviously, we should have J±= 0. This inconsistency arises due to a lack in the mod- 
elization relative to the singular case of infinite wall. Note that, if dealing with a finite 
barrier, (7/2^) |Ar (g + g') = (27,/;i) [Af Q (g + g') / , where g and Q are small. In 
the case of an infinite well, we are in a situation where K tends to infinity. Because of this 
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inconsistency (the infinite well cannot meet the criteria used in our approximations), this 
term should certainly be discarded. The problem can also be circumvented when building 
the function 

$ = T (0 T) = T +0* i= 2 sin (^"^z^ [lAe-'^'^' ] + {lAf e'^'^' [] (B32) 

which properly describes a solution with a spin lying in the plane perpendicular to the DP 
field and for which J = 0. 

APPENDIX C: [100]-ORIENTED BARRIER, ZEROTH-ORDER WAVE- 
FUNCTION COEFFICIENTS 

The continuity of the wave function defined by Eqs. 14.121 and of its derivative at z = 
and 2; = a for the two spin channels provides the following linear system: 



51+^2 + ^2 ^ ^2+ ^2 
1 -r 2 T 2 2K' 2K' 


= Ar 


(Cla) 


K ^ ^ ^ 2K ^ 2K ^ 


- 4-. 


(Clb) 


v42e-^" + 526-^'^ A2e-^'''+ 1326^'" A^e''^'' 
I T 2 IK' 2K' 


= 


(Clc) 


-A2e-^- + B2e^'^ + ^^^e"^'^ + ^^^e^'^ - i^A.e^'^^ 


= 


(Cld) 


-Bi + —A2 ^^2 +^2 + ^2 

2K 2K 


= il 


(Cle) 


Bi A2 B2 A2 + B2 
K' 2K' 2K' 


K' ^ 


(Clf) 


^^26-^'^ - 7^526^'^ + i2e-^'" + S2e^'" - ise'^'^ 
zK zK 


= 


(Clg) 


2K' 2K' K' 


= 


(Clh) 



The coefficients Ai and Ai which define the intensity of the two spin components of the 
incident wave are known (initial conditions). It could be verified that the determinant of 
this system is non zero. We can calculate the eight coefficients Bi, Bi, A2, B2, A2, B2, A3, 
and ^3 from the eight relations, Eqs. ICll We begin to solve these eight equations to the 
zeroth order in C,/Ko or, in other words, by writing C,/Kq = 0. We note the eight equations 
are then divided into two sets: The first four equations are uncoupled to the last four ones. 

The first four equations are related to the spin [1 1]* and write as: 
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Ai = -B1 + A2 + B2 (C2a) 

iqAi = iqBi - KA2 + KB2 (C2b) 

A^e''^" = A2e-^" + 52e^" (C2c) 

ig/lge^^" = -KA2e~^^ + KB2e^'' (C2d) 

and the last four ones are related to the spin [1 — 1]*. The equations are the same by altering 
(Ai, 5i, A2, B2, K) into (Ii, 5i, I2, B2, K'Y 

This is the usual formulation of the tunnel effect. Because Eqs. ICll are written to the 
first order in C,/Ko, we are looking for a solution to the same order. 

To give an example, we look for the results when the incident wave has a spin [1 1]* 
(^Ai 7^ 0, Ai = 0^ Considering Eqs. 12.251 we note that the approximation given by the last 
term of each equation is almost valid as soon as Ka > 2. In Ref. Il3|, K is of the order of 
0.1 which gives a of the order of 20 A in order that the inequality holds, a value which 
is quite reasonable. 

As Ai = 0, this shows that to the zeroth order in ^/Kq, the results may be summarized 

by 





AO) 
— J2 


A2 


= 




_ AO) 
— /s 


A3 


= 


Bi/Ai 




Bi 


= 


52M1 




B2 


= 



(C3) 



where = fj^^ {q, K) and gf"* = gj^^ {q, K) correspond to the standard case (Subsec. IIIPp 



and can be deduced from Eqs. 12.251 This means that, up to the first order in ^/Kq, the 
results are of the shape: 



(C4) 



f^'\q,K,K'),and 

gf'^ = gf'^ {q, K, K') may be equal to zero. In fact, a calculation up to the first order in E,/Kq 



A2/A, 


_ AO) 




I2M1 = i^/K) f^^ 
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where the factors of /j^"* 
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= gfiq,K,K'), ff^ = 
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via Eqs. ICll involves terms of {C,/K) A2 type, which are of second-order in ^/Kq. Therefore 
/W = = and Eqs. [Cl write as 



A,/A, = jf + {^/Kf A,/A, = {i/K) 



(C5) 



B,/A, = gf' + {i/KY gY' B,/A, = {i/K) g\' 

B,/A, = gf^ + {i/Kf gf^ B,/A, = {i/K) gi'^ 

where /f ^ = /f ^ (g, K, K') and gf = gf (q, K, K'). 

Of course if = and Ai 7^ 0, the results are to be inverted, /j^^ ^resp. g^^^^ is 
comparable to, or smaller than, ^resp. gf^^- In Sec. IIV Cl it is shown that /j^^^i is of 
the order of Aj and gf^Ai is of the order of Bj. 
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FIGURES 



FIG. 1. Sketch of the tunnel geometry with definition of notations. The spin-orbit- 
split barrier material of thickness a (medium II) is located between two free-electron like 
materials (media / and ///). The tunnel axis, normal to the barrier, is the z axis. In the 
free-electron-like materials, the real electron wave vector in the z direction is referred to 
as q. In the barrier material, the evanescent wave vector along the z axis is referred to as 
Q + iK, where Q and K are real quantities. The transverse wave vector component, in the 
barrier plane, ^ is conserved in the tunnel process. Then, the overall wave vectors in the 
three media are respectively k/ = k/// = C + Q and k// = ^ + Q+zK. 

FIG. 2. This figure illustrates transformations which, starting from a state of wave vector 
k and with a mean value of the Pauli operator (ct), construct degenerate states. Kq is the 
complex conjugation and K = —iayKo is the Kramers time-reversal operator. K yields a 
state with the wave vector — k* and with the mean spin — (a). The state of wave vector k* 
may be associated to another spin state, corresponding to the mean value (ct)'. Applying 
K to this state, we form a degenerate state with the wave vector — k and associated to the 
mean spin — (ct)'. Four states are finally obtained. 

FIG. 3. Plot of the real-energy lines inside the gap for k = [C,,0,iK], where K and ^ 
are real and positive and tan 6 = ^/K. The calculation is performed using a 14 x 14 k ■ p 
Hamiltonian. The loops are drawn versus ||k|| in 27r/ao units, where is the cubic lattice 
parameter. In all these directions, the spin degeneracy is lifted. Their shape and extension 
sharply depend on 6. For 6 = 43.2°, the two branches are too close to be resolved at this scale. 
The parameters used in the calculation are: P = 9.88 eV.A, P' = 0.41 eV.A, Eq = 1.519 
eV, Ac = Ersc-Enc = 0.171 eV, Px = 8.68 eV.A, A = 0.341 eV, Ea = Enc-Ersc = 2.969 
eV, A' = —0.17 eV (see Ref. |30|for a complete discussion). Inset: Real-band structure (left, 
dashed line. For clarity, only a valence band which is connected to the evanescent branch is 
drawn) and evanescent band across the band gap (right, full line) along the [001] direction 
{6 = 0) where the DP "exchange" field is zero (no spin splitting). 

FIG. 4. Mathematical plot of the real-energy lines for k along [110] as a function of 
the real part of the wave vector Q in the barrier. The calculation is performed for a ratio 
7/7^ = 0.438 A. We are only concerned with negative energies, which refer to evanescent 
states. More precisely, the physical states are located within a very small energy domain 
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below the origin. The domain Q > refers to up-spin states, whereas the domain Q < 
refers to down-spin states. In each case, the imaginary component of the wave vector can 
take the values ±i^. Thus, at a given energy, we have exactly the four possible states 
{Q ± iK) t and {—Q ± iK) |. The down-spin states are Kramers conjugates of the up-spin 
states. 

FIG. 5. This figure is a special case of Fig. 2, when the DP field x li^s along a real 
direction n. Following the same procedure, four degenerate states are constructed, which 
now have their spins quantized along the same direction n (i.e. {a) = (ct)'). 

FIG. 6. (Color online). The lower part of this figure illustrates the spin-dependent tun- 
neling scheme in the case of a [001] oriented barrier (Perel's case). The horizontal plane 
describes the electron wave vector in the barrier; K is taken along the [001] axis and $, lies 
in the barrier plane, along [100]. The upper part of the figure {E > 0) corresponds to the 
real conduction band - the wave vectors are real quantities - and the parabola-like curves 
describing spin-split states along the [101]-direction are drawn. An up-spin state (full line, 
open circle) with the wave vector q' is degenerate with a down-spin state at the wave vector 
q (dotted line, dark circle) and also with up- and down-spin states at the wave vectors — q 
and — q' respectively. This is useful for the calculation of a quantum well, given in Appendix 
IB 2[ Concerning the evanescent states, in a naive effective-mass picture, one may think of 
evanescent states being mirrors of these real states (in the E < domain) with imaginary 
wave vectors. Then up- and down-spin electrons at the energy E, would tunnel with the 
two different wave vectors iq' and iq, thus resulting in a spin-filter effect. However, our 
calculation shows that, concerning evanescent states (lower part of the figure, E < 0), the 
situation is not so simple. In the negative-energy region, the K axis refers to the imaginary 
wave- vector component and $, refers to the real wave- vector component. Real-energy lines 
are found only when tan 6' = ^/K < 1 . These real-energy lines, when drawn for a given 9, 
consist of loops connecting opposite spin states at the zone center ( "up" spin: Full curve; 
"down" spin: Dotted curve. Obviously, when going off the zone center, the spin no longer 
remains a good quantum number - in fact, it can be calculated that its average value rotates 
along the loop - but it has to be pointed out that, in the D'yakonov and Perel' description, 
the energy eigenvectors are also pure spin states which depend on the 9 ratio). Two of these 
loops are drawn here. Let us consider a tunneling process at the energy E (horizontal grey 
plane or yellow plane in the online edition) of an electron with the wave vector component 
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£ in the barrier plane, which has to be conserved in the tunnehng process. It can be ob- 
served here that the two states marked on the loops by a dark circle (K') and an open circle 
(K) - which are energy degenerate - are associated to the same real wave-vector component 
^. However, they correspond to two different 9 as they are respectively associated to the 
imaginary components iK and iK', along the tunneling direction. The difference between 
K and K' results in a spin-filter effect. Inset (upper left): Top view of the plane at energy 
E showing the intercepts with the loops which determine the relevant wave vectors K and 
K'. 
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